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At stronger gauge-field couplings, the domain wall fermion (DWF) residual mass (wres), a mea- 
sure of chiral symmetry breaking (Xsb)^ grows rapidly as it is largely due to near zero fermion 
ly-j eigenmodes of log(r), where T is the 4D transfer matrix along the fifth dimension. The number 

CNJ of these eigenmodes increase rapidly at strong coupling. To suppress such near-zero eigenmodes, 

^vq we have added to the DWF path integral a multiplicative weighting factor consisting of a ratio of 

determinants of Wilson-Dirac fermions having a chirally twisted mass with a large negative real 
component and a small imaginary chiral component. Numerical results show that this weighting 
factor with an appropriate choice of twisted masses significantly suppresses mies while allowing 
adequate topological tunneling. 
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It is known that the unwanted near-zero eigenmodes of the Wilson Dirac operator can be sup- 
pressed by including a suitable multiplicative weighting factor in the lattice path integral. Vranas [1, 
2] demonstrated that a determinant of Wilson- Dirac fermions with real mass equal to the domain 
wall height acts as such a weighting factor for domain wall fermions (DWF). Fukaya [3] demon- 
strated that a ratio of determinants Wilson Dirac fermions with a real numerator mass and a chirally- 
twisted denominator mass acts as such a weighting factor for overlap fermions. 

Here, we have selected as a weighting factor for DWF a ratio of determinants of Wilson Dirac 
fermions where both numerator and denominator have chirally-twisted masses. The first section 
reviews properties of mres useful for interpreting the action of the weighting factor; the second 
section characterizes this weighting factor, and the third section presents numerical studies of this 
weighting factor. Readers interested in our numerical results can skip to the last section. 

Note that this report goes beyond the previous conference presentation. At the time of the 
conference, numerical studies has just begun and the reported results were of limited scope. At this 
time, numerical studies have progressed and more interesting results can be reported. 

1. Relevant Properties of DWF 

XsB in domain wall fermions (DWF) can be represented in an effective low-energy Lagrangian 
as an extra fermion mass term with mass wires, where m^s can be defined as a ratio of the matrix 
element of the axial current at the 5D mid-plane and the matrix element of pseudo-scalar density. 
A useful approximation to mres is the following [4]. 



Here, p(A) is the density of eigenmodes of a 4D Hamiltonian, H/[ = — log(r), where T , a 4D 
transfer matrix acting along the 5'th dimension, is a functional of the Wilson-Dirac operator with 
negative mass, D^{—M) [5]. The term exponential in arises from eigenmodes with eigenvalues 
I A I greater than Ac (mobility edge) and have extended 4D support and exponential decay in the 5'th 
dimension. This term is easily controlled by increasing Lg. 

At larger and stronger gauge-field couplings, the second term proportional to L^' (the 
inverse term) is larger and of principal interest. This term arises from eigenmodes with near-zero 
eigenvalues (in any case | A | < A^) having localized 4D support. Note that p (0) is used to represent 
the density of near-zero eigenmodes, not of eigenmodes exactly at 0. 

It is believed that many (perhaps all) 4D-localized, near-zero eigenmodes arise from patches 
of enhanced roughness in the gauge fields that are localized to a few lattice spacings. Such patches 
have been associated artifactually with the appearance or disappearance of instantons in Monte 
Carlo time. Gauge-field roughness, and therefore p(0) and mres, increases rapidly as couplings, j8, 
becomes stronger, and increase particularly rapidly at couplings below the QCD transition. 

There is of course a caveat: Proper sampling of gauge-field topology requires that the MD 
trajectories tunnel between topological sectors, and such tunneling is always associated with sign- 
change of a Dirac eigenmode. MD trajectories therefore must not be prevented from passing 
through gauge fields where one or more Dirac eigenvalues are very close to (exactly in the 
continuum). The very-close-to-zero eigenmodes relevant to topological tunneling are also thought 
to be 4D-localized. 
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Figure 1: Behavior of the 20 lowest eigenvalues of Dp-{—mo) versus mo at a lattice spacing offl-^ «2.0Gev 
for, (a), a pure Wilson gauge action, and, (b), a pure Wilson gauge action combined with two Wilson-Dirac 
fermions. 

2. The Auxiliary Determinant 

Arguably the ideal weighting factor to suppress m^ss would be built from H4, which is compu- 
tationally inaccessible. Thus, our and others' weighting factors are built from the computationally- 
inexpensive Hermitian Wilson-Dirac operator Y=,Dy^{—M) M > 0. As the determinant of this 
operator is near zero if one or more eigenmodes are near zero, including the squared determinant 
as a multiplicative factor in the DWF path integral will under-weight corresponding regions of 
gauge-field phase space during Monte Carlo sampling, or equivalently, generate a repulsive force 
away from such regions during molecular dynamics (MD) evolutions. 

We point out that this use of the Hermitian Wilson-Dirac operator relies on the assumption: 
the eigenmodes operator y=,D'^{—M) are 'sufficiently similar' to the those of the operator 7/4 so 
that suppressing the near-zero eigenmodes of the former operator also suppresses the near-zero 
eigenmodes of the latter operator. Similarity is to be expected, because, in the case of a continuous, 
unbounded 5'th dimension. Ha, become the Hermitian Wilson-Dirac operator, and, when this two 
operator has a zero eigenmode, the Hermitian Wilson-Dirac operator does also. Vranas has pro- 
vided a proof of principle that using a determinant of Wilson-Dirac fermions suppresses near-zero 
eigenmodes [2]. Figure 1 from that paper demonstrates this suppression. 

However, the Wilson fermion determinant is not an ideal weighting factor: it essentially pre- 
vents topological tunneling since it is very close to zero in the presence of very-close-to-zero eigen- 
modes; and it over-weights regions of gauge-field phase space where all eigenmodes are large since 
it is large in such regions. 

Our weighting factor avoids these problems by including a small imaginary addition to the 
mass in the numerator (see refs. [6] and [2]) and the denominator, the latter proposed in ref. [3]. 



W{M,ef,eb) 



det[a^{-M + iehfyDw{-M + iebY^)] 
det[D^{-M + iefY^YD^{-M + ieff)] 
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Figure 2: versus |A| for gy/e,, of 0.001/0.10, 0.005/0.50, and 0.040/0.50. Exemplary val- 
ues at |A,| of 0.001, 0.025, and 0.100 of ^,(0.01,0. 10)/^,(0.04, 0.50) are: 16, 2.9 and 0.6; and of 
jri(0.005,0.50)/^,(0.04,0.50) ai-e 62, 3.4 and 1.2, respectively 



^ dti[Dy^{-MyD^{-M)]+e} Xl + e} 
det[D^(-M)tDy (-M)] + el ^1 xf + el 

Where A, are the eigenvalues J'Dyf. The effect of can be more readily appreciated by con- 
sidering the total force, ^ , and the force from a single eigenmode, generated during an MD 
simulation. ' 

^te..)^-^(-.„*O^Ej^(-.n*0|^(|^)_~E^.te,..)i(|^)^, (2.3, 

T is the MD time; = (A? + ef) /(A? + ; and usually < < £^ < 1. 

We examine in detail the variation of with assuming that {dXi/dU) varies sufficiently 
slowly with for the three combinations of Ef/Eb for which numerical results are presented 
below: 0.001/0.10, 0.005/0.50, and 0.040/0.50. Figure 2 illustrates the variation of with 
\Xi\ and Usts exemplai-y values of the ratios ^,(0.01,0.10)/^,(0.04,0.50) and ^,(0.005, 0.50)/ 
^,(0.04,0.50) for of 0.001, 0.025, and 0.100. This figure makes apparent that the size and 
position of the maximum of the force and the size and size and position of the tail of the force can 
be significantly and independently controlled by the imaginary mass parameters Ef and £/,. Table 1 
summarizes the | A,- [-dependence of relative forces and eigenmode suppression. 

3. Numerical Studies 

The numerical results reported here have had the goal of rapidly obtaining a qualitative under- 
standing of the properties of our selected weighting factor, and therefore have been continued only 
for sufficient trajectories, generally from 300 to 500 (MD timestep usually 0.2), to obtain, first, a 
reasonably accurate values for mres^ second, an estimate the QCD transition temperature, Tc (15c), 

'This expression for ^ abstracts features useful for the subsequent discussion. A more correct derivation could 
proceed according to an ansatz introduced in [7]: first differentiate the MD Hamiltonian with respect to MD time, and 
then find the force as the traceless, anti-Hermitian (TA) part of the coefficient of the MD momentum, is a formally 
color-independent coefficient and factors out of the TA operation. 
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Table 1: Dependence of relative forces and eigenmode suppression on \Xi\. 

For 0.1 < the forces satisfy ^,(0.005,0.5) « ^,(0.04,0.5) > 2 * ^,(0.01,0.1); and 
accordingly the weighting factors ^(0.005,0.5) and ^(0.005,0.5) sup- 
press eigenmodes in this range equally and to a greater extent than the 
weighting factor ^(0.01,0.1). 
For < 0.025 the forces satisfy ^,(0.005,0.5) > 2 * ^,(0.01,0.1) > 4 * ^,(0.04,0.5); 

and accordingly the weighting factors suppress eigenmodes in this range 
inthe order /T (0.005, 0.5) > #'(0.01,0.1) > ^(0.04,0.5). 



and third, an approximation of the topological charge. Reference data without the weighting factor 
were drawn from prior studies of domain wall thermodynamics. All simulations use the Iwasaki 
gauge action, have a weighting with two Wilson fermions in the numerator and two in the denomi- 
nator, have two light DWF (masses: 0.03), and one heavy DWF (mass: 0.037) for = 32. Further 
simulations details are the same as in [4] (except that this reference used = 16). 

Physical scales have not, to date, been measured, nor has the susceptibility of Accord- 
ingly, we compare individual data values from different simulations only at the estimated jSc, and 
compare only groups of values above and below jSc. Further, j8c has been estimated by visual char- 
acteristics of the fluctuations of versus trajectory that were found to be reliable proxies for the 
susceptibility in prior thermodynamic studies. ^ We interpret the data in Tables 2 and 3 in view of 
Figure 2 and the discussion in Section 1. The data are consistent with the following hypothesis. 

njres suppression Depends on suppression of near-zero eigenmodes of the Hermitian 

Wilson-Dirac operator having |A | values less than about 0.1 — 0.5; 
Topological charge sup- Depends on suppression of very-close-to-zero eigenmodes of the 
pression Hermitian Wilson-Dirac operator having | A | values less than about 

0.01-0.05; 

First, compare the expected suppression of eigenmodes with 0.1 < | from Table 1 with the mpes 
measurements from Table 2. It is expected that, in comparison to mres without any weighting factor, 
ffires is most Suppressed at 0.005/0.50; suppressed to a somewhat less degree at 0.04/0.50; and 
clearly less suppressed at 0.01/0.10. As expected, measured mjes, in comparison to mies without 
any weighting factor, is most suppressed at 0.005/0.50 and less suppressed at 0.01/0.10. mres 
suppression at 0.04/0.50 is intermediate between that at 0.005/0.50 and at 0.01 /O.IO. ^ 

Further, comparing mres without any weighting factor to mjes at 0.01/0.10 above and below 
the transition: below the transition, relative suppression of mres is greater than at the transition; and 

^Two such characteristics are relied on. First, in the vicinity of j8e, the amplitude of 'J"!', fluctuations is a maximum. 
Second, the nature of these fluctuations differs above and below the transition. For j3 > j3c, 'PI' has upward, short, 
relatively large-amplitude, fluctuations (spikes) without corresponding downward fluctuations superimposed on a largely 
constant baseline. Spike amplitude and frequency increases as fic is approached. For /3 < /3e, the baseline disappears and 
fluctuations become substantially symmetrical, and include low amplitude fluctuations extending over 10-50 trajectories. 

"'/Ores is known to be largely independent of temperature, as is confirmed by comparing mj-es at Ef/Cb = 0.005/0.50 
and /3 = 1.75 from the zero temperature lattice 16** x 32 with that from the non-zero temperature lattice 16"' x 8 x 32. 
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Table 2: Values of Wres for various Ef, Sb, and j3. On the left are thermodynamic runs for a series of j3 at 
different combinations of and £/,. The horizontal double line through this sub-table identifies the estimated 
transition, and divides the table into an upper part where j3 < jic and a lower part where j3 > jSc- On the right 
are zero-temperature results for the series of £/ at a selected j3 and e^. (* denote values that were determined 
from fewer than 300 trajectories.) 



lattice 16^ X 8 X 32 


no weighting factor 


£f/£b -- 


= 0.01/0.10 




= 0.005/0.50 


P 






"^res 




"^res 






1.875 


0.0101(5)* 










1.900 


0.0072(3)* 






1.95 


0.0252(1) 


1.925 


0.0054(4)* 






2.00 


0.0102(1) 


1.950 


0.0035(3)* 


1.750 


0.0015(1) 


2.05 


0.0046(2) 


1.975 


0.0026(4)* 


1.800 


0.0009(3)* 


2.08 


0.0022(2) 


2.000 


0.0020(2)* 


1.850 


0.0003(4)* 


2.11 


0.0011(1) 










2.14 


0.0007(1) 











lattice 16"^ X 32 


j3 = 1.75 £fc = 0.50 






0.040 


0.0025(1)* 


0.020 


0.0018(1) 


0.005 


0.0014(10)* 



Table 3: Values of the average of the absolute topological charge, < \ Vfop\ >, for various Ef, Et,, and j3. The 
left and right sub-tables have the same structure as in Table 1. Runs were too short to determine useful error 
bounds; it is believed that errors where < \Vtop\ >— 0{l) may be up to 50%, and where < \Viop\ 0(0.1) 
may be up to 100%. (* denote values that were determined from fewer than 300 trajectories.) 



lattice 16^^ X 8 X 32 


no determinant 


£f/£b = 


= 0.01/0.10 




= 0.005/0.50 




< \Vtop\ > 


i8 


< \v,op\ > 




< \ Vtop\ > 






1.875 


1.6* 










1.900 


1.7* 






1.950 


2.7 


1.925 


1.0* 






2.125 


1.7 


1.950 


0.8* 


1.75 


0.1 


2.050 


0.9 


1.975 


0.5* 


1.80 


0.2* 






2.000 


0.4* 


1.85 


0.0* 



lattice 16"^ X 32 


= 1.75 Eh = 0.50 




< \ Vtop\ > 


0.040 


2.1* 


0.020 


0.1 


0.005 


0.4* 



above the transition, relative suppression is less than at the transition. Accordingly, this weighting 
factor can be expected to be increasingly effective as gauge-field coupling increases. 

Next, compare the expected suppression of eigenmodes with | A, | < 0.025 from Table 1 with the 
topological charge measurements from Table 3. It is expected that, in comparison to topological 
charge without any weighting factor, topological charge is most suppressed at 0.005/0.50; less 
suppressed at 0.01/0.10; and minimally suppressed at 0.04/0.50. Again as expected, measured 
topological charge, in comparison to topological charge without any weighting factor, is virtually 
eliminated at 0.005/0.50 and suppressed to an intermediate degree at 0.01/0.10. Although data is 
limited, the topological charge at 0.04/0.50 is likely to be closer to the topological charge without 
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any weighting factor than to the topological charge at 0.01/0.10. 

We conclude that near-zero eigenmodes of the Hermitian Wilson Dirac operator can usefully 
serve as proxies for the near-zero eigenmodes of H^, and also that the near-zero and very-close-to- 
zero eigenvalue ranges are reasonable for the characteristics of the DWF used here. It is likely that 
the eigenvalue ranges that here appear relevant to suppression of wires and of the topological charge 
are specific to the parameters defining the DWF fermions simulated here, particularly the explicit 
light quark masses. For example, in the case of a different explicit light quark mass, the eigenvalue 
ranges relevant to suppression of the topological charge can be expected to change"*. However, 
the eigenvalue ranges relevant to suppression of m^es can be expected to be little changed, as 7/4 is 
independent of explicit quark mass. 

Nevertheless, the parameters studied here may well be close to those needed to allow the 
use of DWF for larger lattice spacing and consequently on larger lattice volumes with /Wres values 
better controlled that has been possible to date. At present the largest physical volumes that have 
been accessible are L 2.7fm boxes used with 24^ x 64 lattices at \/a = 1.73Gev and 32^ x 64 
lattices at 1/a = 2.42Gev work of RBC/UKQCD with m^s values of 0.00315(2) and 0.000676(11), 
respectively. Pending a more accurate determination of the lattice scale, the results presented here 
suggest that DWF simulations with mres ~ 0.002 and 1/a = 1.4fm giving a 32^ lattice volume 
with a linear size of L = 4.5fm may be possible. Such larger physical lattice volumes will offer 
important opportunities for study of QCD thermodynamics for study of ^ ^ ;r7r decays. 
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the transfer matrix, is independent of mf, and det[£)£ny/r(l)] is the Pauli-Villars subtraction. 
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Appendix 

Subsequent to submitting this report to the Proceedings of Science, we determined the behavior 
of several lowest Wilson-Dirac eigenvalues versus mass for sample configurations from the runs 
reported in Table 2, and present in this Appendix results for thermalized sample configurations 
from the three runs appearing above the double line in the left portion of Table 2.^. These three 
runs are expected to be below, but near to, the QCD transition. When no weighting factor is used, 
the transition has been independently determined to be at or near jS = 2.031(5) ^. 

These results are sufficient to qualitatively appreciate the effects of the weighting factor. How- 
ever, the statistics are insufficient for quantitative analysis at this time. In all cases, the domain wall 
fermion height was -1.8. 
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Figure 3: Without weighting factor: behavior of the 12 lowest eigenvalues of D-0^{mo) versus mo at for a 
thermalized 16^^ x 8 x 32 lattice at j3 = 1.950. 
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(b) j3 = 1.750 and £f/£i, = 0.005/0.50. 



Figure 4: With weighting factor: behavior of the 12 lowest eigenvalues of D^{mo) versus niQ at for a 
thermalized 16^ x 8 x 32 lattice at ^3 = 1.950 and Ef/eb = 0.01 /O.IO (left subfigure) and for j3 = 1.750 and 
Ef/eh = 0.005/0.50 (right subfigure). 



^Diagrams of this type are often presented — otq > on the horizontal axis. 

^Michael Cheng et al., The transition temperature using 2+1 flavors of domain wall fermions at Nt = 8, in prepara- 
tion. 
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